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I. INTRODUCTION 



Quantum entanglement, as a kind of distinctive quantum correlation without the classical 
counterpart, is widely recognized as the crucial resource of quantum communication and 
computation hi. However, there are also other nonclassical correlations that can even exist 
in separated quantum states [21, I3|. The finding that nonclassical correlations other than 
quantum entanglement may provide the speedup in the deterministic quantum computation 
with one pure qubit (DQCl) protocol j^, |5[ has greatly motivated the related study. 

One of the essential issues is to distinguish different kinds of correlations in quantum 
systems. In the classical information theory, correlation is perfectly characterized by the 
Shannon entropy, which is represented by the classical mutual information 6|. When it 
comes to quantum systems, the Shannon entropy is replaced by the von Neumann entropy 
and the total correlation in a bipartite system is characterized by the quantum mutual 
information [7] . For bipartite systems, the quantumness of correlation can be quantified by 
the quantum discord which represents the difference between classical information theory 
and quantum information theory 2]. A related method concerning classical correlation is 
proposed and it is based on the maximal information one can extract with a one-sided 
local measurement [3|. It has been demonstrated that almost all quantum states contains 
non-vanishing quantum discord that is composed by the quantum entanglement and non- 
entanglement quantum correlation [8|. 

There are many investigations of quantum correlation measured by quantum discord 
in different kinds of physical systems, including the spin chains 9l-|l2|. Jaynes-Cummings 
systems 13|, spin-boson systems [l4l], optical systems 15|, quantum dots [16] and NMR 
systems 



systems 
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Recently, this kind of measurement is also extended to continuous variable 



- |20| and even the quantum biology systems ^J. Quantum discord has been 
found to be useful in quantum information processing and quantum information theory. 
The non-entanglement quantum discord may provide the quantum advantage in the DQCl 



quantum computation protocol j^, ISj. The vanishing of quantum discord implies completely 



positive dynamics maps |22|. Quantum discord also plays important roles in many basic 
physical problems, for example the Maxwell's demon |23,], quantum phase transition J9l-|l2| 



and the relative effect 
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25| 



The quantification of quantum discord is based on the one-sided measurement on a bi- 



partite system, which is usually asymmetric. We would get different results by choosing 
different subsystem to be measured. A symmetrical method with two-sided measurement 
over both partitions of a bipartite system is proposed to quantify the classical correlation 
which is represented by the maximal classical mutual information 26|, |27|]. There are also 
other approaches proposed to distinguish classical and quantum correlations. Inspired by 
the work consumption when extracting information from a heat bath, a thermodynamics 



approach is used to defined quantum correlation in quantum systems |28|, |29|. In partic- 
ularly, the difference between the extractable information with the closed local operation 
and classical communication and the total information is defined as the quantum informa- 
tion deficit, which is used to quantify quantum correlation |29[. Because classical states are 
measured without disturbance, classical and quantum correlations are also characterized by 
the measurement-induced disturbance [Sy]. Recently, by employing the relative entropy as 
a distance measure of correlations, the method proposed by Modi et al. [31] provides an 
unified view of quantum and classical correlations. Different from quantum discord based on 
bipartite mutual information, this method can be extended to quantified different kinds of 
correlations in multipartite systems of arbitrary dimensions. Generally, all these measure- 
ments of classical and quantum correlations mentioned above are not equal to each other 
due to the different definitions. However, they are consistent in the case of distinguishing 



classical and quantum correlations in Bell diagonal states 15|. 

Another interesting subject is to investigate the dynamics of different kinds of correlations 
in noisy environments. For one side, the inevitable interaction between a quantum system 
and its environment would destroy correlations in the system and leads to the reduction of 
useful resource. The knowledge of the dynamic behavior of correlations will help us to design 
suitable protocols to protect correlations under processing. For the other side, stimulated 



by the discovery of distinctive dynamic behavior of entanglement sudden death [321 , that is, 
entanglement disappear completely in a finite evolution time, the investigation of the unusual 
dynamic behavior of classical and quantum correlations has caused great interests. Usually, 
the dissipative correlation evolution is essentially dependent on the types of noises that act on 
the system. Markovian noises would cause the irreversible decay of system information into 
the environment. It has been shown that quantum discord is more robust than entanglement 



under Markovian noises 33|, in the sense that entanglement may suffer from sudden death. 



whereas quantum correlation decays exponentially. The sudden change in behavior in the 



decay rates of classical and quantum correlations has been predicted in different kinds of 



Vlarkovian environments 3J] and has been experimentally verified in an optical system 

15|. The decoherence-free evolution of quantum and classical correlations under certain 

Markovian noise and the sudden transition from classical to quantum decoherence regime 

are also observed 
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35| . The Markovian dynamics of classical and quantum correlations 
is also experimentally investigated in a nuclear magnetic resonance quadrupolar system 
[l7|. When it is extended to the non-Markovian environment with memory effect, the 
feedback information from the environment to the interested system may greatly affect the 
dynamic behavior of correlations. Quantum discord under non-Markovian noises is shown 
to be instantaneously vanished compared to the completely disappearance of entanglement 
in a finite time interval |36l . |37|. The revival of quantum discord without the revival of 



entanglement is also demonstrated in non-Markovian noises 10 1. 

In this paper, based on the equal footing method with the relative entropy as a distance 
measurement of correlations, we experimentally investigate the dynamics of classical and 
quantum correlations of a two-photon Bell diagonal state in a one-sided non-Markovian 
environment, which is simulated by a Fabry- Perot cavity followed by quartz plates. At the 
beginning^ evolution, the sudden transition from classical to quantum decoherence regime 
is shown 
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35| . Due to the refocusing effect of the relative phase in the non-Markovian 
environment, the quantum correlation revives from a near zero area and then decays again 
in the subsequent evolution. However, the non-Markovian effect is too weak to revive the 
classical correlation, which remains constant at the same evolution range. We then perform 
a ax operation on the photon under decoherence and the sudden transition from quantum 
to classical revival regime is obtained, in which correlation echoes are formed. 

The paper is organized as follows. The measure of correlations using relative entropy and 
the theoretical description of the experiment are given in section II. The experimental setup 
and results are presented in section III. We then give a discussion and conclusion in section 
IV. 



II. RELATIVE ENTROPY OF CORRELATIONS AND THEORETICAL DE- 
SCRIPTION OF THE EXPERIMENT 

A. Relative Entropy of Correlations 

The magnitude of a specified property in a quantum system can be quantified by the 



distance from the interested state to the closest state without that desired property |31 |. 
For example, quantum entanglement can be characterized by the relative entropy of en- 
tanglement {REE) [38], which is described as the minimal distance measured with relative 
entropy between the state p and a separated state a and is expressed as REE = minS'fpl Icr). 
S{p\\a) = — tr(plog2cr) — S{p) and S{p) = — tr(plog2p) = —J2jPj^'^S2Pj i^ the von Neu- 
mann entropy {pj represent the eigenvalues of p). D is the set of separable states. As a 
result, if all the distances are measured with relative entropy, different kinds of correlations 
in a quantum system can be measured on an equal footing. The quantum correlation is 
then defined as the minimal distance between p and a classical state x? which is expressed 



as Q = minS'(pllx) [31]. C represents the set of classical states. Whereas the classical corre- 
lation is defined as the minimal distance between x and a product state vr and is expressed 
as C = minS'(x||7r) [31], where P is the set of product states. 

TTgP 

It has been demonstrated that the calculations of quantum and classical correlations in 



a quantum system (p) can be further simplified as [31 1 



Q = Sixp) ~ S{p), (1) 

C = S{n.J-S{x,), (2) 

where Xp represents the closest classical state of p and n^^^ is the corresponding reduced state 
of Xp in the product form which is the closest product state of Xp- 
The total mutual information of p can be calculated as 

X = 5(7r,)-5(p), (3) 

which represents the minimal distance between p and its reduced state in the product form 



TTp [311] • For bipartite systems, X is equal to the quantum mutual information. 

Generally, it is difficult to find the closest classical state Xp- However, for the Bell diagonal 
state p = J2i=i ^i\^i){^i\ where Aj are the non- increasing eigenvalues and |$i) represent the 



four Bell states, the analytic expression of Xp is found and Xp = ^^'^'^^ (|$i)($i| + |$2)($2|) + 
^^^(|$3)($3| + |$4)($4|) [31]. As a result, for the Bell diagonal state with the eigenvalues 
{Ai, A2, A3, A4}, the eigenvalues of the closest classical state is {|(Ai + A2), |(Ai + A2), ^{Xs + 
A4), |(A3 + A4)}. We can therefore calculate the quantum correlation Q according to equation 
([1]). The product states of p and Xp are both identical and equal to the normalized identity 
1/4. As a result, the classical correlation and total mutual information can be calculated 
according to equations (^ and (^ respectively and X = Q + C. The analytic solution of the 



relative entropy of entanglement for Bell diagonal states is given by [38 1 

REE = 1 + Ai log2 Ai + (1 - Ai) log2(l - Ai), (4) 

if Ai > 1/2, whereas 

REE = 0, (5) 

if Ai G [0,1/2]. 

As a result, with the knowledge of p, one can compute different kinds of correlations. 

B. Theoretical description of the experiment 

The pure dephasing environment is a kind of uniquely quantum noise, which causes 
randomness between the relative phases of information carries. In optical systems, the 
coupling between the photon polarization states (information carriers) and photon frequency 
(noise freedoms) in a birefringent environment leads to the dephasing with a trace over 



frequency [39 1. 



Consider a maximally entangled polarization state of two photons with the form of |\1') = 
^{\HH)a^b + \HV)a,b + \VH)a^b — \W) a,b) , whcrc \H) and \V) represent horizontal and 
vertical polarization states, respectively. The subscripts a and b denote different paths of 
these two photons. When the entangled state passes through dephasing environments, which 
are simulated by quartz plates with the optic axes set to be horizontal, the final polarization 
state for a certain single frequency can be written as 

1$') = ^{\HH)a,b + e''^''\HV)a,b + e'*^\VH)a,b 

-e*^'^"+'^'')|rV),,b), (6) 



where 0a = La^nua/c and 0b = Lb/S.nujh/c. La (uja) and L^ {ujh) represent the thick- 
ness of quartz plates (the frequency of the photon) in paths a and 6, respectively. An is 
the difference between the indices of refraction of ordinary and extraordinary light, and 
c represents the vacuum velocity of the photon. By tracing over all the frequency de- 
grees of freedom, two decoherence parameters k^ = j g{uja)e^'^''dua and i^b = f /(wb)e*'^''dcJh 
would impose on the reduced polarization density matrix, where g{u}a) and f{uJb) repre- 
sent the frequency distributions of the photon in paths a and b and they are normalized as 
f g{ua)dijJa = 1 and f f{ub)dijJb = 1, respectively. The final density matrix in the canonical 
basis {\HH), \HV), \VH), \VV)} becomes 



m 
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P=A 



1 

Kn 



K„ 



K Kh 



^a^b * 



f^al^b 



1 -Kt 



\ l^af^b f^a f^b 



(7) 



/ 



where k* (k^) corresponds to the complex conjugate of Kq (k;,). This final evolved state can 
be transformed into a Bell diagonal form with local unitary operations. For a special case 
with Kq and Kb both setting to be real, the four eigenvalues of p are given by {^(1 + /ta)(l + 
Kfe), \{l-Ka){l + Ki), \{l + Ka){l-Ki,), i(l -K^) (1 -Kfe)}. Becausc < «„ < 1 and < k^ < 1, 
the maximal eigenvalue is |(1 + /t(j)(l + Kb) and the minimal eigenvalue is |(1 — /ta)(l — ^^b)■ 
Whereas the identification of the second maximal eigenvalue is dependent on the relative 
magnitudes of Kq and k^. If we set /€„ to be a fixed value (with a fixed thickness of La) 
representing the decoherence parameter in preparing the initial mixed state and k^ to be the 
decoherence parameter in the evolution ranging from 1 to 0, the second maximal eigenvalue 
is ^(1 — Ka){l + Kfe) when Ka < Kf,, whereas the second maximal eigenvalue is given by 
^(1 + Ka){l — Kb) when Ka > Kb- Therefore, the four non- increasing eigenvalues of the closest 
classical state Xp are {^(1 + Kb), ^{l + Kb), |(1 — Kb), \{l — Kb)} when Ka < Kb. In the case 
of Ka > Kh, the four eigenvalues become {^(1 + Ka), \{1 + Ka), |(1 — Ko), ^(1 — Ka)}, which 
are all fixed values. As a result, the quantum correlation is calculated as 



Q 



i(l + Ka) log2(l + Ka) + Ul- Ka) log2(l 



Kn 



if Ka < Kb, 



Ml + Kb) log2(l + Kb) + 5(1 - Kb) log2(l - Kb), if Ka > Kb- 



(8) 



and the classical correlation is expressed by 

/- _ J li^ + l^b) ^Og2{l + Kb) + 1(1- l^b) ^Og^il- Kb), if Ka<Kb, 
1^1(1 + Ka) log2(l + Ka) + |(1 - ^^a) loggll " /^a), if ^a > H- 

We can find that the quantum correlation and classical correlation remain constant when 
Ka < Kb and Ka > Kb, respectively (for fixed Ka). They overlap at the point of Ka = Kb. As 
a result, the state of equation ([7]) with Kb = 1 represents the kind of initial states with the 
property of exhibiting the sudden transition from classical to quantum decoherence regime 

mhm- 

Generally, the frequency spectrum of the photon is peaked at a central value Uq with 
a finite width a, for example the Gaussian function like frequency distribution f{ub) = 
{2/ ^/^^a) exp[—4:{uJb — ujo)^/o-'^]. The decoherence parameter Kb is therefore calculated as 
Kb = exp[— {Lb An/ c)'^a'^/ 16 + i{LbAn/c)uo] and it decays exponentially, which leads to 
the Markovian limited dynamics of correlations. However, if the frequency distribution 
of the photon in mode b becomes discrete, such as the combination of finite A^ Gaussian 
frequency distributions /(wt) = J2j=i^j{'^/V^'^j)^w[~^{^b — ^j)'^/'^]] where Aj are the 
relative amplitude for each Gaussian function distribution with the central frequencies Uj 
and frequency widths aj. In this case, the decoherence parameter is calculated as Kb = 
^ ■ ^ y4j exp[— (LfeAn/c)^cr|/16 + i{LbAn/c)uj]. During the dephasing process, the overall 
relative phase may refocus and the non-Markovian effect occurs, which leads to the revival 
of Kb. Kb may be larger than Ka again and it would give rise to the revival of correlations. In 
experiment, the discrete frequency distribution can be realized by passing the photon in mode 
b through a Fabry- Perot cavity, which behaves as an optical resonator [4:U\. Wavelengths 
for which the cavity optical thickness is equal to an integer multiple of half wavelengths are 
resonant in the cavity and transmitted. Other wavelengths within the reflective band of the 
Fabry-Perot cavity are reflected. 

By controlling the non-Markovian effect, we can get the revival of Q without the revival 
of C (the maximal revival value of Kb is less than Ka), and the case that both Q and C get 
revival (the maximal revival value of Kb larger than Ka). Actually, the maximal revival of Q 
and C can be realized by completely compensating the decoherence in the pure dephasing 
environment. If we exchange the polarizations of \H) and \V) of the photon in mode b 
during the dynamics, the randomness of the relative phase caused in the previous evolution 



time is compensated by the same subsequent evolution time [39|. For the state |\&'), it 
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becomes |^') = le'''"'{\HV)a,b + \HH)a,b + e''^''\VV)a,b - e''''''\VH)a,b)- As a result, the effect 
of decofierence in mode b reduces to an unobservable global phase and the final state is 
changed to the initial form with the exchanging of \H) and \V) in mode b again. In our 



experiment described below, we employed such spin echo like technology 4l| to obtained 
the maximal revival of both Q and C. 

III. EXPERIMENTAL SETUP AND RESULTS 

Figure [1] shows the experimental setup to investigate the correlation dynamics in the 
non-Markovian environment. The second harmonic ultraviolet (UV) pulses are frequency 
doubled from a mode- locked Ti: sapphire laser with the center wavelength mode locked to 
0.78 /im (with a 130 fs pulse width and a 76 MHz repetition rate). These UV pulses are 
then focused into two thin, identically cut type-I /3-barium borate (BBO) crystals with 



their optic axes aligned perpendicularly to each other |42|]. Degenerate photon pairs with 
wavelengthes centred at 0.78 yum, created from the spontaneous parametric down conversion 
process, are emitted into paths a and b. By compensating the birefringence effect in the 
BBO crystals with quartz plates (CP) in both paths, the prepared maximally entangled 



state {l/V2{\HH)a,b + \VV)a,b)) can has a high purity |43|. 

A half-wave plate (HWPl) with the optic axis set at 22.5°, which changes \H) into 
1/V2\H + V) and \V) into l/\/2\H — V), is used to transfer the maximally entangled state 
into the exact form of |\I'). Quartz plates (Ql) with the optic axis set to be horizontal are 
used to dephase the photon in the path a. The frequency distribution g{uJa) is considered 
as a continuous Gaussian function, which is defined by the interference filter with a 3 nm 
full width at half maximum (FWHM). 

The frequency distribution of the photon in path b {f{ub)) becomes discrete after it passes 
through the Fabry-Perot (FP) cavity, which is a 0.2 mm thick quartz glass with coating films 



(reflectivity 90% at wavelengths around 780 nm) on both sides 40|. Quartz plates (Q2 and 
Q3) with the optic axes set to be horizontal are used to introduce the dephasing effect. A 
half-wave plate (HWP2) with the optic axis set to 45° acting as a a^ operation can exchange 
\H) and \V). Another half- wave plate (HWP3) with the same setting as HWP2 after Q3 is 
used to change the final form of the output state. 

The density matrix of the final state is reconstructed by the quantum state tomography 



process [44], where the 16 different coincidence measurement bases are set by quarter- wave 
plates (QWP), half- wave plates (HWP) and polarization beam sphtters (PBS). Both photons 
are then detected by single-photon detectors equipped with 3 nm interference filters to give 
coincidence counts. 
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FIG. 1: (Color online). Experimental setup. Maximally entangled polarization photon pairs, 
created from the spontaneous parametric down conversion process by pumping the two-crystal 
geometry type I BBO with ultraviolet (UV) pulses, are emitted into paths a and b. These two 
photons pass through quartz plates (CP) to compensate the birefringence in BBO. A half-wave 
plate (HWPl) and quartz plates (Ql) in path a are used to prepared the required mixed state. The 
Fabry-Perot (FP) cavity together with quartz plates (Q2 and Q3) and two half-wave plates (HWP2 
and HWP3) in path b is employed to monitor the dephasing coupling. After passing through 
quarter- wave plates (QWP), half- wave plates (HWP) and polarization beam splitters (PBS) which 
allow to tomographically reconstruct the density matrix, both photons are then registered by 
single-photon detectors (Dl and D2) equipped with 3 nm interference filters (IF). 

In our experiment, the final output state is a Bell diagonal state with the form of equation 
dlj). According to the equal footing method [SlJ, experimental results of quantum correlation 
(Q), classical correlation (C) and total correlation (X) are deduced from equations ([1]), ([2]) 
and (|3]) with S{7i^^) = ^(vrp), respectively. The relative entropy of entanglement (REE) 
can be calculated from the maximal eigenvalue of p according to equations (jl]) and ([5]) |38 |. 
As a result, all kinds of correlations can be deduced from the density matrix of p. 

Fig. |2] displays the dynamics of correlations, as a function of the thickness of Q2 which 
is represented by the corresponding retardation x {L = x/An). The thickness of Ql rep- 
resented by the retardation is set to be llTAo (Ao=0.78 /im is the central wavelength of 
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FIG. 2: (Color online), (a) Experimental results of the non-Markovian dynamics of correlations. 

Black squares, dark cyan diamonds, blue stars and red dots represent the experimental results 

of X, C, Q and REE with the black solid line, dark cyan dashed line, blue solid line and red 

dotted line representing the corresponding theoretical prediction, respectively. The inset (the x 

axes represent the total thickness of Q2) represents the corresponding evolution of eigenvalues of 

p and Xp- The magenta, dark cyan, blue and dark yellow dots represent the experimental results 

of Ai, A2, A3 and A4 with the magenta, dark cyan, blue and dark yellow solid lines representing 

the corresponding theoretical predictions, respectively. Black squares and red stars represent the 

experimental results of ^(Ai + A2) and |(A3 + A4) with the black and red solid lines representing the 

corresponding theoretical predictions, respectively, (b) Theoretical predictions of the correlation 

dynamics with the FWHM of f{oJb) identically fitting to 0.2 nm. The black solid line, dark cyan 

dashed line, blue solid line and red dotted line represent X, C, Q and REE, respectively. Aq = 0.78 

/im. 



the photon) to prepare the initial mixed state and k^ is equal to about 0.607. The HWP2 
and HWP3 are not used {i. e., without a^. operation) and Q3 is set to 0, in which the 
phenomenon of entanglement collapse and revival occurs [40(] . We find that there exhibits 



a sudden transition from classical to quantum decoherence area 



35( 1 ■ At the beginning 



11 



of the evolution, the quantum correlation Q (blue stars) remains constant and then decays 
exponentially after the thickness of about 120Ao. Due to the refocusing effect of the relative 
phase, Q revives from near zero at the thickness of about 440Ao and reaches the maximum 
of 0.11 at about 540Ao. With further increasing Q2, Q decays monotonically again. The 
classical correlation C (dark cyan diamonds) behaviors quite differently. It decays exponen- 
tially at the beginning and then remains constant all the time after the thickness of 120Ao 
despite the non-Markovian effect. Q and C overlap at the thickness of 120Ao, in which 
the sudden change in behavior in their decay rates are observed 15|, |3J]. The evolution of 



relative entropy of entanglement REE (red dots) is also shown, which suffers from sudden 



death J32|] at the thickness of about 189Ao and is consistent with our previous results 
If we continue to increase Q2, REE also revives to its maximally value at the thickness of 
about 540Ao (the value is relative small in the figure). The evolution of total correlation 
X (black squares) first decays exponentially and then revives just as that of Q. The black 
solid line, dark cyan dashed line, blue solid line and red dotted line represent the theoret- 
ical predictions of X, C, Q and REE. The inset displays the corresponding dynamics of 
eigenvalues of p and Xp- The magenta, dark cyan, blue and dark yellow dots represent the 
experimental results of Ai, A2, A3 and A4 with the magenta, dark cyan, blue and dark yellow 
solid lines representing the corresponding theoretical predictions, respectively. Whereas the 
black squares and red stars represent the experimental results of |(Ai + A2) and |(A3 + A4) 
with the black and red solid lines representing the corresponding theoretical predictions. It 
can be seen clearly that the sudden transition of classical and quantum decoherence occurs 



at the point when the switch in the second maximal eigenvalue A2 occurs 35[ and it is 
consistent with our previous theoretical prediction. At the period with the relative phase 
refocusing, the four eigenvalues Aj behave correspondingly, i. e., Ai (A3) increases and A2 
(A4) decreases. However, |(Ai + A2) and |(A3 + A4) remain constant and there is not revival 
of C according to its definition. The errors of the experimental results are estimated using 



the method proposed in ref. 4J], which is mainly due to the random fluctuation of each 
measured coincidence counts (the errors from the uncertainties in aligning the wave plates 
is relatively small). In this approach, the errors are directly deduced from the fluctuation 
of the corresponding eigenvalues. The error bars of the relative entropy of entanglement 
involve only the maximal eigenvalue of p, which are relative small in the figure. 

In our experiment, the frequency distribution f{ujb) is treated as three Gauss-like wave 
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packets centered at 778.853, 780.160 and 781.459 nm with the relative probabihties of 0.37, 



0.44 and 0.19, respectively 40|. The FWHM of these wave packets are identically considered 
as 0.85 nm and we obtain good fittings. The maximal revival value of k^ is about 0.385, 
which is smaller than /«„. As a result, the classical correlation remains constant after the 
sudden transition point in the case of Fig. [2^, which is immune from the non-Markovian 
effect and is consistent with our previous analysis in the theoretical part. With the narrower 
FWHM of f{uJb), the maximal revival value of Kb can be larger than Ka and the switch in 
the second maximal eigenvalue would occur again, in which the refocusing effect would be 
strong enough to revive the classical correlation. Fig. |2]d shows the theoretical prediction of 
the correlation dynamics with the FWHM of /(wf,) identically fitting to 0.2 nm (the maximal 
revival value of kj, is about 0.944). We can see that X (black solid line), C (dark cyan dashed 
line), Q ( blue solid line) and REE (red dotted line) are all revived. The sudden transition 
from quantum to classical revival regime is obtained at the thickness of about 477Ao (the 
revival value of Kb is equal to Ka). 

It is a great challenge to realize the small FWHM to obtain the observable revival of 
classical correlation in the experiment. We then implement a a^ operation on the photon 
in path b to investigate the correlation dynamics and get the revival of classical correlation. 
The thickness of Q2 is first increased to 200Ao or 400Ao followed by a cTx operation which 
is fulfilled by the HWP2. The thickness of Q3 is then increased from zero without further 
increase of Q2 to get the corresponding dynamics of correlations, which are shown in Fig. |3k 
(Q2=200Ao) and Fig. [3b (Q2=400Ao), respectively (the x axes represent the total thickness 
of Q2 and Q3). The black squares, dark cyan diamonds, blue stars and red dots represent 
the experimental results of X, C, Q and REE with the black solid line, dark cyan dashed line, 
blue solid line and red dotted line representing the corresponding theoretical predictions, 
respectively. We can see that correlation echoes are formed when Q3=200Ao in Fig. [3^ and 
Q3=400Ao in Fig. |3b, i- e., at the time when Q3 is increased to the same thickness as 
Q2. This implies that the dephasing effect between \H) and \V) caused by Q2 is completely 
compensated by Q3 and all the correlations maximally revives to the initial values. During 
this process, the sudden transition from quantum to classical revival area is observed (Fig. |3^ 
at the thickness of about 280 Aq and Fig. [3]d at the thickness of about 680 Aq). The classical 
correlation revives from the constant period to its maximal value, whereas the quantum 
correlation stops revival and remains constant. If we further increase Q3, the subsequent 
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FIG. 3: (Color online). Experimental results of correlation echoes: (a) with the ax operation at 

200 Aq; (b) with the ax operation at 400 Aq- Black squares, dark cyan diamonds, blue stars and 

red dots represent the experimental results of X, C, Q and REE with the black solid line, dark 

cyan dashed line, blue solid line and red dotted line representing the corresponding theoretical 

predictions, respectively. The x axes represent the total thickness of Q2 and Q3. The blue arrows 

identify the ax operation points. The insets in (a) and (b) (the x axes represent the total thickness 

of Q2 and Q3) represent the corresponding evolution of eigenvalues of p and Xpi respectively. The 

magenta, dark cyan, blue and dark yellow dots represent the experimental results of Ai, A2, A3 and 

A4 with the magenta, dark cyan, blue and dark yellow solid lines representing the corresponding 

theoretical predictions, respectively. Black squares and red stars represent the experimental results 

of 2('^i + ^^2) and 2 ('^3 + -^4) with the black and red solid lines representing the corresponding 

theoretical predictions, respectively. Aq = 0.78 ^m. 

dynamics of correlations are just the same as that shown in Fig. [2^. The blue arrows in 
the figures identify the ax operation points. The correlation dynamics between the initial 
state and the maximally revived states are symmetric about the ax operation point, which 
is similar to the phenomenon of spin echo in nuclear magnetic resonance 4l(| . The insets in 
Fig. |3^ and Fig. [3b represent the corresponding dynamics of eigenvalues of p and xp- The 
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magenta, dark cyan, blue and dark yellow dots represent the experimental results of Ai, A2, 
A3 and A4 with the magenta, dark cyan, blue and dark yellow solid lines representing the 
corresponding theoretical predictions, respectively. Whereas the black squares and red stars 
represent the experimental results of |(Ai + A2) and ^(As + A4) with the black and red solid 
lines representing the corresponding theoretical predictions, respectively. We can see that 
the eigenvalues also display the echo effect and the second switch in the second maximal 
eigenvalue A2 in both the insets of Fig. |3^ and |3]d leads to the revival of classical correlation. 

IV. DISCUSSION AND CONCLUSION 

The reasons for quantum advantage in quantum information processing are still contro- 



versial 45|. Recently, a simple one-to-one relationship between bipartite entanglement of 



formation [46| and quantum discord in a general tripartite system is proposed 47| . By tak- 
ing the environment which is initially maximally entangled to the system with the reduced 
maximally mixed state in the DQCl protocol into consideration, it is suggested that both 
ihe quantum discord and entanglement are responsible for the quantum computer speedup 



48 



47 1 . Inspired by the factorization law of entanglement evolution in noisy quantum channels 



50|, it is expected that similar simple relationship for quantum correlation under noise 
environment exists. With the increasing interests, there will be more distinctive discoveries 
in this field. 

In our experiment, we have investigated the correlation dynamics in a non-Markovian 
dephasing environment. The whole point of a non-Markovian environment is that it retains 
a memory of a system at a given time and then later passes this information back into the 
system in some form or another. In this experiment, the non-Markovian environment acts via 
the FP cavity followed by quartz plates on the biphoton system of only one of the photons. 
Due to the discrete frequency distribution of the photon which leads to the refocusing of 
the relative phase in the dephasing environment, the non-Markovian effect occurs and the 
revival of quantum correlation is obtained from near zero area. However, the non-Markovian 
effect is too weak to revive the classical correlation. With the narrower FWHM of the wave 
packets of the discrete frequency distribution, the revival of classical correlation would be 
also achieved. On the other hand, if the FWHM of the wave packets become lager, the 
non-Markovian effect becomes weaker. When the discrete frequency distribution tends to 
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continuous Gaussian distribution, we would get the Markovian limited dynamics. We further 
implement a cTx operation on the photon in path h to investigate the corresponding correlation 
dynamics. During this process, we obtain correlation echoes, in which the sudden transition 
from quantum to classical revival regime is observed. This work is a useful and informative 



add-on to our previous works, in which the entanglement collapse and revival is observed [40 1 
and the Markovian limited dynamics of correlations is demonstrated [l5|. The method can 
be used to control the revival time of correlations, which would find important applications 
in quantum memory. 

We thank Cheng-Hao Shi for help in performing the experiment. This work was sup- 
ported by National Fundamental Research Program, National Natural Science Foundation 
of China (Grant Nos. 60921091, 10734060, 10874162, 11004185), China Postdoctoral Science 
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